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Correlating anomalies of the microwave sky: The Good, the Evil and the Axis 
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At the largest angular scales the presence of a number of unexpected features has been confirmed 
by the latest measurements of the cosmic microwave background (CMB). Among them are the 
anomalous alignment of the quadrupole and octopole with each other as well as the stubborn lack of 
angular correlation on scales > 60°. We search for correlations between these two phenomena and 
demonstrate their absence. A Monte Carlo likelihood analysis confirms previous studies and shows 
that the joint likelihood of both anomalies is incompatible with the best-fit A Cold Dark Matter 
model at > 99.95% C.L. Extending also to higher multipoles, a common special direction (the 
'Axis of Evil') has been identified. In the seek for an explanation of the anomalies, several studies 
invoke effects that exhibit an axial symmetry. We find that this interpretation of the 'Axis of Evil' 
is inconsistent with three-year data from the Wilkinson Microwave Anisotropy Probe (WMAP). 
The data require a preferred plane, whereupon the axis is just the normal direction. Rotational 
symmetry within that plane is ruled out at high confidence. 
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m ■ 
> ■ 

\o . 

(N ■ 

m ■ 

o : 
r» . 
o ■ 

^ ' 

Or 
i 

o ■ 
a: 



I. STANDPOINT 

With the emergence of more and more precise and de- 
tailed cosmological observations the inflationary A cold 
dark matter (ACDM) model remains to provide a sur- 
prisingly good fit. Thereby the most precise and dis- 
tinguished lever arm is provided by the cosmic mi- 
crowave background (CMB). The standard inflationary 
model predicts approximately scale- invariant, statisti- 
cally isotropic and gaussian temperature fluctuations on 
the surface of last scattering and is fully consistent with 
the data. But after the release of three years of mission 
data from the Wilkinson Microwave Anisotropy Probe 
(WMAP) 0,0 there remain at least open questions and 
at most serious challenges upon the inflationary ACDM 
model of cosmology. 

Based on the high precision measurements of WMAP, 
a couple of anomalies on the microwave sky have been 
identified. These anomalies manifest themselves at the 
largest angular scales, mainly among the quadrupole and 
octopole (the dipolc is overwhelmingly dominated by our 
local motion with respect to the CMB), but also extend- 
ing to somewhat higher multipoles. The corresponding 
anomalies may be divided into two types: 

First, and already seen by the Cosmic Background 
Explorer's Differential Microwave Radiometer (COBE- 
DMR) 0] and confirmed by the first-year analysis of the 
WMAP team 0, there is a lack of angular two-point cor- 
relation on scales between 60° and 170° in all wavebands. 
In Q the angular two-point correlation function of the 
three-year WMAP measurements has been computed. 
Going form COBE-DMR to WMAP(3yr) the lack of cor- 
relation persists and moreover it has been outlined [5| 
that among the two-point angular correlation functions 
none of the almost vanishing cut-sky wavebands matches 
the reconstructed full sky and neither one of the latter 
matches the prediction of the best-fit ACDM model. This 
disagreement has been shown to be even more distinc- 



tive in the WMAP(3yr) data than in the WMAP(lyr) 
data and is found to be unexpected at 99% C.L. with 
respect to the three-year Internal Linear Combination 
[ILC(3yr)] cut-sky. Recently, it has been shown Q that 
indeed quadrupole and octopole are responsible for the 
lack of correlation and that most of the large-scale an- 
gular power comes from two distinct regions within the 
galactic plane (only 9% of the sky). 

Second, there exist anomalies concerning the phase re- 
lationships of the quadrupole and octopole. There are a 
number of remarkable alignment anomalies 0, @], e.g. an 
unexpected alignment of the quadrupole and octopole 
with the dipole and with the equinox at 99.7% C.L. and 
99.8% C.L., respectively @. In contrast to such extrinsic 
alignments, that is alignments of the low multipoles with 
some physical direction or plane, like the dipole or the 
ecliptic, the intrinsic alignment between quadrupole and 
octopole does not know about external directions. In this 
work, we address the intrinsic alignment of quadrupole 
and octopole with each other, which from the ILC(3yr) 
map is found to be anomalous at 99.6% C.L. with re- 
spect to the expectation for an statistically isotropic and 
gaussian sky |5|. 

Both types of CMB phenomena challenge the state- 
ment of statistical isotropy of the CMB sky at largest 
angular scales. Here we want to study the relation be- 
tween the lack of angular correlation and intrinsic align- 
ment of quadrupole and octopole. 

In Q it has been shown that intrinsic alignments 
among multipole moments extend also to higher ones and 
it has been proposed that the strange alignments at large 
angular scales involve a preferred direction, called the 
'Axis of Evil'. This axis points approximately towards 
(l,b) ~ (—100°, 60°) and is identified as the direction 
where several low multipoles (£ = 2 — 5) are dominated 
by one m-mode when the multipole frame is rotated into 
the direction of the axis. Recently, the analysis of 
the 'Axis of Evil' has been redone in the light of the 
WMAP(3yr) with the use of Bayesian techniques [ll|. It 
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was argued [T3| that the 'Axis of Evil' is rather robust 
against foreground contaminations and galactic cuts. A 
recent [l3| cross-correlation of CMB data and galaxy sur- 
vey data shows no evidence for an 'Axis of Evil' in the 
observed large scale structure. In contrast, recently an 
opposite claim has been put forward [T3 |. 

Motivated by these observed CMB anomalies, several 
mechanisms based on some axisymmetric effect have been 
proposed, although the operational definition of the 'Axis 
of Evil' [3, [l(| does not necessarily imply the existence of 
such a strong symmetry. Among the various effects that 
have been suggested to possibly introduce a preferred 
axis into cosmology are: a spontaneous breaking of sta- 
tistical isotropy [151 ] , parity violation in general relati vity 
(l6| . anisotropic perturbations of dark energy lal. 
residual large-scale anisotropies after inflation BJ l20| . 
or a primordial preferred direction (2l| . At the same 
time it has been studied [22|, [H, [HI how the local Rees- 
Sciama effect [25j of an extended foreground, non-linear 
in density contrast, affects the low multipolc moments 
of the CMB via its time-varying gravitational potential. 
In a scenario with a single overdensity the coefficients of 
the spherical harmonic decomposition, the ag m , become 
modified by only zonal harmonics, i.e. m — modes. 
This is equivalent to an axial effect along the line con- 
necting our position with the centre of the source. 

In fact, the observed pattern in the CMB for 
quadrupole and octopole is a nearly pure au mode re- 
spectively; as seen in a frame where the z-axis equals 
the normal of the plane defined by the two quadrupole 
multipole vectors (26|. In 0] it has been argued already, 
that foreground mechanisms originating from a relatively 
small patch of the sky would mainly excite zonal modes. 
Moreover all additive effects where extra contributions 
are added on top of the primordial fluctuations would 
have difficulties explaining the low multipole power at 
large scales without a chance cancellation. 

In this work we study how the inclusion of a preferred 
axis compares with the intrinsic multipole anomalies at 
largest scales. Our analysis is restricted to axisymmetric 
effects on top of the primordial fluctuations from stan- 
dard inflation, thus secondary or systematic effects. We 
quantify how poorly an axisymmetric effect at low mul- 
tipoles of whatever origin matches the three year-data 
of WMAP. We demonstrate that there is no correlation 
between the two types of intrinsic \ow-£ anomalies: the 
two-point correlation deficit and intrinsic alignment; and 
that there remains none even when a preferred axis is 
introduced to the problem. 

The paper is organized as follows. In the next sec- 
tion the necessary technical framework is introduced, es- 
pecially the multipole vector formalism. Thereafter in 
Sec. IIIII we recapitulate the best- fit standard inflationary 
model ('the Good') and its quantitative predictions for 
the \ow-£ microwave sky as well as the anomalous ('the 
Evil') measurements of WMAP(3yr). Before we conclude 
in Sec. El the consequences of an axial symmetry ('the 
Axis') on the low-£ CMB sky are discussed in Sec. IIVI 



II. CHOICE OF STATISTIC 

A common observable is the multipole power. Accord- 
ing to the standard perception of inflationary cosmology, 
the CMB fluctuations are believed to follow a gaussian 
statistic and to be distributed in a statistically isotropic 
way. The notion of statistical isotropy means that the ex- 
pectation value of pairs of coefficients {o^, m ia^ m ) is pro- 
portional to 5e.'t.5m'm- Usually the proportionality con- 
stant measuring the expectation value of the multipole 
on the full sky is estimated by 
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with Ti being the £-tYi multipole of the CMB temperature 
anisotropy. It can be expanded with the help of spheri- 
cal harmonics as: Tg = Y] m ag m Yi m . Note that, since we 
consider multipole moments that are real, the ag m must 
fulfill an additional condition: a* tm — (— l)™ l a^_ m . Us- 
ing the estimator ([1} the angular two-point correlation 
function is given by 
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where the Pg are the Legendre Polynomials of l-th order. 

Besides of the multipole power itself, it is useful to in- 
troduce an all-sky quantity that embraces all scales. As 
inspired by the S1/2 statistic, presented in Q for measur- 
ing the lack of angular power at scales larger than 60°, 
we use here an analogous all-sky statistic [5j: 
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This is a measure of the total power squared on the full- 
sky. In contrast to the Sx/2 statistic [J] the S{ u ii statistic 
does not contain any a priori knowledge on the variation 
of the two point angular correlation for angles > 60° . 
Here we are considering especially the large angular scales 
but we are not interested in the monopole and dipolc and 
thus arrive at 
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Of course all multipoles have to be considered for the 
full-sky statistic §5§ but we use the truncated part ([J|, 
because here the anomalies are most pronounced and we 
want to check for the interplay of this part of the full-sky 
power statistic with the other (phase) anomalies within 
quadrupole and octopole. This part is then simply to 
be added to the rest of the sum of (squared) multipole 
power in ([3]), recovering the expression for the full-sky. 

Next we turn to the statistics involving the phase 
relationships of multipoles. We use the concept of 
Maxwell's multipole vectors [27| in order to probe sta- 
tistical isotropy, since this representation proved to be 
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FIG. 1: Mollweide projection of the sky with quadrupole (upper row) and octopole (lower row) multipole vectors [Equation 
J5J]. The mesh consists of steps in 30°. Displayed are ten pairs of quadrupole vectors (small dots) and their ten area vectors 
[Equation ([6| (big dots)] as well as ten triples of octopole vectors (small dots) and their area vectors (big dots); togetherness 
is indicated by color. The arbitrary sign of the vectors has been used to gauge them all to the northern hemisphere. The 
statistically isotropic and gaussian case (left column) is broken by the imprint of a strong axial effect a^o = 1000/iK (right 
column) whereupon multipole vectors move to the pole and area vectors move to the equatorial plane. The onset of the shown 
separation of multipole vectors and cross products can already be observed at moderate axial contributions of ato ~ 100/iK. 



useful for analyses of geometric alignments and special 
directions on the CMB sky. Normally the CMB data is 
decomposed into spherical harmonics and the coefficients 
dim containing the physics. Alternatively, with the use 
of the multipole vectors formalism [28| we can expand 
any real temperature multipole function on a sphere into 

i 

Tt(6,<p) = ^ ai m Y hn (8,ip) 



= A® 



(5) 

Therein the z-th vector belonging to the multi- 

pole I and e(6,(p) — (sin 8 cost/?, sin 8 sin <p, cos 8) stands 
for a radial unit vector. With this decomposition all of 
the information in a multipole moment is reorganized 
such that we obtain a unique factorisation into a scalar 
AW counting the total power and £ unit vectors v^ 1 ^ 
encoding all the directional information. The product 
expansion term alone contains also terms with 'angular 
momentum' ^ — 2,^ — 4, ... These residuals are subtracted 
with the help of the term Ce(8, <p). Because all the signs 
of the multipole vectors may be absorbed in the scalar 
quantity A^> the multipole vectors are unique up to their 
sign which carries no physical meaning. For definiteness, 
we define them to point to the northern hemisphere. 



In order to disclose correlations among the multipole 
vectors we first consider for each I the £(£ — l)/2 in- 
dependent oriented areas built from the cross products 
SI, HI: 
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whereof we will also use the normalized vectors n^' 1 '^ = 
w (t;hi) /\ w (Z;i,j)^ Now, in [1] and subsequent works the 
dot products of the area vectors proved to be a handy 
expression in order to quantify alignments of the mul- 
tipole vectors among each other and also with external 
directions (which we do not consider here). The follow- 
ing measure as stated in 29] and used in 5, 8, 26] serves 
as a natural choice of a statistic in order to quantify the 
intrinsic alignment of quadrupole and octopole oriented 
areas: 
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Note that we consider only the very largest scales, i.e. we 
use the statistic only for £ = 2, 3. Analogously, a statistic 
involving the normalized area vectors is given by: 
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FIG. 2: Evolution of the Monte Carlo likelihood of the alignment statistics S nn © an d S ww ([8}. The effect of an axis in 
the CMB is modeled via increasing additional zonal harmonics with coefficients aeo- At a^o = 1000/iK the multipoles become 
purely zonal in good approximation. Regarding WMAP's ILC(3yr) map S n n is unexpected at 98.3%C.L. and S ww is odd at 
99.5%C.L. with respect to the statistically isotropic and gaussian sky (bold histograms). The best improvement is reached for 
both statistics at roughly aeo = 100/iK. 



III. INFLATIONARY ACDM PREDICTIONS 

Standard inflationary ACDM cosmology requires the 
CMB anisotropics to be gaussian and statistically 
isotropic. For the subsequent analysis we have produced 
Monte Carlo realisations of the harmonic coefficients aim 
following the underlying ACDM theory. From [28] an al- 
gorithm is available which we use to obtain Monte Carlo 
multipole vectors from the coefficients. Mollweide maps 
of a sample of random gaussian and statistically isotropic 
quadrupole and octopole vectors as well as their normals 
are given in Figure [1] (left column) . 

Concerning the question of correlations between the 
multipole power and the alignment of multipole vectors, 
it appears natural to expect that there is none. That 
is because we invoked gaussian random and statistically 
isotropic skies leading to multipole vectors ([5|) indepen- 
dent of the multipole power ([T]) . This assumption needs 
to be tested and quantified. 

Nevertheless, a small correlation could be expected 
from the following reason: Considering only multipoles 
up to some limiting power, the resulting probability den- 
sity distribution for the a^ m must be non-gaussian. In 
fact, this restriction leads to a negative kurtosis for the 
ai m distribution (the skewness vanishes). Having that in 
mind, it appears suddenly unclear whether the naive ex- 
pectation of vanishing correlation of power with intrinsic 
alignment will hold. Below we substantiate the absence 
of correlations by means of a Monte Carlo analysis. 

Let us first look at the alignment anomalies. In Figure 
[2] the likelihood of the quadrupole and octopole align- 
ment statistics S ww and S nn is shown. The predictions 
of the standard inflationary ACDM model are shown as 
the bold histograms respectively (= vanishing axial con- 
tamination). According to the three-year ILC map from 
WMAP |2| we get the following measured values for the 



alignment statistics: 

5 ILC(3yr) = _ 86g2 > = 0.7604 

when [f| corrected for the Doppler-quadrupole. The to- 
tal number of Monte Carlos we produced per sample 
is N — 10 5 . We infer that the unmodified inflationary 
ACDM prediction is unexpected at 98.3% C.L. with the 
Snn statistic and unexpected at 99.5% C.L. [30] with 
respect to the S ww statistic. 

Next we consider the cross-correlation between the 
intrinsic phase anomalies and the multipole power ([T]) 
within the low-£. For this we chose those ag m that al- 
low for say the lowest possible 5% in the left tail of the 
distributions for C2 and C3 that follow from statisti- 
cal isotropy, gaussianity and the ACDM best-fit to the 
WMAP data. Then we compute the expression S ww 
for the selected ai m and compare it to the according 
ILC(3yr) value. As expected, no correlation is found, 
that is neither the shape nor the expectation value of 
the alignment statistic is shifted. We find the same also 
for the combination of the lowest allowed 5% in C2 and 
the highest 5% from the right tail of the distribution of 
C3 and the remaining two possible combinations thereof. 
As we do not find any correlations, we can conclude that 
the S ww and S nn statistics are not sensitive to the non- 
gaussianity induced by the restriction to low multipole 
power. 

Moreover, we probe the opposite direction by tagging 
those ai m that lie in the allowed right tail of the S ww dis- 
tribution with respect to The distribution of 
the multipole power for C2 and C3 made of these ai m re- 
mains unchanged. The latter finding confirms that mul- 
tipole power and the shape of multipoles (phases) are 
uncorrelated. 

Using Equation (|U), the [2| Maximum Likelihood Es- 
timate (MLE) from the WMAP ILC(3yr) map for the 
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sample size N 


joint p 


error A 


100000 


0.048% 


0.008% 


100000" 


0.001% 


0.002% 


10000 


0.02% 


0.02% 


1000 


0% 


0.02% 



"With respect to WMAP(lyr) ILC data 



TABLE I: Joint likelihoods © for S^T" and S ww being in 
accordance with data simultaneously. The experimental val- 
ues refer to WMAP's ILC(3yr) map [2j except for the second 
row. The error A of the factorisation in equation @ is the 
difference between left hand side and right hand side in that 
equation. 

angular power spectrum yields S^^ c ' = 29431/iK 4 . 
Compared to the value of 136670>K 4 from the ACDM 
best- fit to WMAP(3yr) data, this is not significantly un- 
expected, with an exclusion level of only 92.1% C.L. 

Now we want to check for correlations between the all- 
sky multipole power and the multipole alignment. As 
for reasons explained in the next section we prefer the 
S ww statistic to S nn in the following correlation analy- 
sis. In Figure [3] the scatter plot of S ww against Sj. ™ c 
is shown. The form of the contour can be understood 
as just the folding of the x 2 -like form of the distribution 
for S"f ™i" c with the gaussian-like form of the S ww distri- 
bution. At first glance we see from Figure [3] that the 
MLE from WMAP(3yr) ^uac.MLE = 2 g43i M K 4 requires 
the alignment statistic to be of middle values (around 
0.4), which is inconsistent with the respective measured 
anomalous value from ILC(3yr). Moreover the lack of 
any linear behaviour in the contour suggests that there 
is no correlation between the two statistics. 

Given that no correlation is present between S ww and 
S'f™ 1 nc , we would expect that the joint probability that 
both power and alignment are in accordance with data 
factorizes according to: 

P (SfuT C < data A S ww > data) = 

Pi (Sf3T < data ) P2 (S ww > data) . (9) 

But in reality we can only access finite statistical sam- 
ples of these quantities and the factorisation will not be 
exact. However we will check the validity of within 
our statistical ensemble. When using the full sample 
with N = 10 5 respectively we obtain a joint likelihood 
of p ~ 0.05%. The error A of the factorisation, which we 
define as the difference between the left hand side in © 
and the right hand side, is of the order 0(1O~ 5 ), that is 
of the order of the Monte Carlo noise. In order to track 
the evolution of the error A we also compute the joint 
likelihood © for smaller subsamples; see Table [TTTJ Re- 
ducing N to N — 10 we obtain an even smaller joint 
likelihood of p = 0.02% but with an error that is of the 
same magnitude. With N = 10 3 we do not have a single 
hit for the joint Monte Carlos leading to p = 0% with 
the same error as in the N — 10 4 case of A = 0.02%. 
Note that just one Monte Carlo hit in favor of the joint 
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FIG. 3: Contour of the scatter of intrinsic alignment (0 ver- 
sus full-sky power squared Q. The shape can be understood 
from the folding of the two respective distributions. The to- 
tal number of Monte Carlo points is N = 10°. The mea- 
sured data point from WMAP three-year data is included. 
The maximum of likelihood requires S W w far smaller than 
obtained from ILC(3yr). Consistency with the data can be 
excluded at 99.95% C.L. Contours correspond to lines of 1/2" 
times the maximal likelihood, with n — 1, ... ,5. 

case would raise the error here to A = 0.08%. In the end, 
the convergence of the joint likelihood appears to be very 
slow with respect to the sample size N. 

Furthermore we are interested in the stability of the 
results for A with respect to changes in the measured 
data. For this we choose the WMAP(lyr) values: 

£trunc,pscudo-C f = ^54^4 &nd ^C(lyr) = Q 773^ 

We use a sample of the full size N = 10 5 and obtain 
a joint likelihood with respect to the one-year data of 
p = 0.001% with an error A = 0.002%. That is, with 
respect to one-year data both the joint likelihood and its 
error are of the order of the Monte Carlo noise. From 
the WMAP(lyr) data alone we could exclude the joint 
case (j9]) rather conservatively at 99.99% C.L. This ap- 
pears to be a stronger exclusion than the one from three- 
year data. But we do not bother much about the differ- 
ence because of the different estimators that have been 
used by the WMAP team for the angular power spectrum 
(pseudo-C £ vs. MLE) 0. 

We quote the most conservative result, namely the full 
sample joint likelihood case for S ww and S?^ji nc with re- 
spect to the WMAP(3yr) data. Therefore we can exclude 
that case at > 99.95% C.L. with an error in the third digit 
after the comma lying within the Monte Carlo error of 
the used sample (N = 10 5 ). 

Finally we analyze the correlation of the all-sky power 
statistic S'fuJJ 110 and the intrinsic multipole alignment S ww 
by quantitative means: 

It is well known from statistics, that when checking a 
finite two-dimensional sample for correlations, the em- 
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FIG. 4: The sign of additional axial contributions apo has no physical effect on the statistics S nn and S W w ■ For the quadrupole 
this follows from the symmetry of the Legendre Polynomial P2 [see Equation \V2f\. The quadrupole contribution is kept fixed 
at d2o = 100/xK while the axial contribution to the octopole is varied both in magnitude and in sign. Respective pairs of 
±<23o histograms lie virtually on each other and their statistics are thus indistinguishable. The reference histograms following 
from the axially unmodified ACDM model (bold histograms in Figure [2} lie nearly on top of the displayed 0,20 ~ lOOtiK and 
d3o = ±10/iK cases, and are thus not shown. 



piric covariance 
cov[ 5fuil nc > S ww } = 
1 N 

_ T (^fMM _ ^furi" ) (S W W,i — S ww ) (10) 

i=l 

is a crucial quantity. The bar stands for the mean of a 
variable. As the covariance is a scale dependent measure, 
i.e. depending on the magnitudes of the sample values 
S W w,i and S ww _i, the dimcnsionless Bravais- Pearson co- 
efficient or empirical correlation coefficient is the better 
expression to use: 

cov[S^ nc , S ww ] 

fuU ' WW V C0V [ ^full ' Sful\ nC ] C0V [ 

(11) 

Finally, employing the WMAP(3yr) data we obtain an 
empirical correlation coefficient of Ps^'^.S^ = —0.0027 
with respect to the full sample N = 10 5 , which indeed 
indicates only marginal correlation. 

IV. INCLUSION OF AN AXIS 

Now we ask what happens when introducing axial con- 
tributions on top of a statistically isotropic and gaussian 
microwave sky. The presence of a preferred direction with 
axisymmetry in the CMB will exclusively excite the zonal 
modes in case the axis is collinear to the 2-axis. Here we 
do not bother about external directions since the internal 
alignments are independent of these. Therefore such an 
axis will manifest itself through additional contributions 
a£Q. We are considering the quadrupole and the octopole 



and the question arises, in how far the sign of the axial 
contributions ±a«) plays a role. The coefficients a£ m can 
be reconstructed from 

f AT 

a im = J —(9,v)Y; m <m. (12) 

Obviously, within the quadrupole the sign of ±0,20 is ir- 
relevant because of the symmetry of the Legendre Poly- 
nomial P2 with respect to 6 = 90°. The Legendre 
Polynomial P3 however is antisymmetric with respect to 
9 = 90°. Therefore the relevance of the sign of the oc- 
topole contributions 030 has to be clarified. Consequently 
we have chosen a fixed value for the axial quadrupole 
contribution a 2 o and have then varied the according oc- 
topole contribution in sign and in magnitude. The results 
are displayed in Figure [4] Apparently the S nn and S ww 
statistics that are important here, do not distinguish be- 
tween the sign of the applied axial effect. Therefore we 
need not to bother about the signs of the a^o and let 
them henceforth be positive. 

In Figurc[5]thc evolution of the S ww and S nn statistics 
with respect to increasing axial contributions is displayed 
in terms of likelihood histograms: 

Let us first look at the evolution of the S nn statistic. 
This expression measures the average | cos | of the angles 
between the quadrupole oriented area and the octopole 
areas. The pure Monte Carlo peaks at 0.5 reflecting the 
fact that the average distance of four isotropically dis- 
tributed vectors on a half-sphere from each other is 60° 
in the case of statistical isotropy. It is a half-sphere be- 
cause the signs of the multipole vectors are arbitrary and 
so we choose them all to point to the northern hemi- 
sphere. When increasing the contribution of the axial 
effect the multipoles become increasingly zonal and ar- 
rive at being purely zonal in a good approximation at 
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FIG. 5: Scatter contour of pairs of S»„ and Sfun" after an 
axial modification of am = 70/iK has been applied; this is 
the contribution involving maximal improvement in S-w-w (see 
Figure^. The total number of Monte Carlo pairs is N = 10 5 . 
Note that the horizontal axis now runs from zero to 1.4 x 
10 _6 mK 4 , whereas in Figure [3] the maximal displayed value 
is 4 x 10~ 7 mK 4 . The inclusion of a preferred axis leaves all-sky 
multipole power and intrinsic alignment totally uncorrelated 
and inconsistent with the WMAP(3yr) data. Contour lines 
are defined as in Figure [3] 



values of ago = lOOO^K. On the level of the multipole 
vectors this means that their cross products all move to 
the equatorial plane (see Figure [TJ . That is the reason 
why the histogram in Figure [2] (left) moves to the right 
when we increase the axial effect, because now isotropy is 
broken from the half-sphere to the half-circle making the 
Snn histogram peak sharper at higher values. The mea- 
sured value from the ILC(3yr) map of ^ L „ c(3yr) =0.868 is 
anomalous at 98.3% C.L. with respect to the pure Monte 
Carlo (bold histogram in Figure [5]) which stands for the 
statistically isotropic and gaussian model. By adding ax- 
ial contribution the maximal improvement is reached at 
am = 100/iK where the ILC(3yr) becomes unexpected 
at 96.7% C.L. Further enhancement of the axial effect 
makes the S nn statistic more and more narrow around 
an expectation value < 0.7. This makes it impossible 
to remove the anomaly in the S nn cross-alignment with 
respect to the ILC(3yr) experimental value only by in- 
creasing the axial contribution to high enough values. 

On the other hand the S ww statistic additionally mea- 
sures the modulus of the sin of the angles between the 
multipole vectors themselves. As can be seen from Fig- 
ure [T] multipole vectors are all moving toward the north 
pole clustering more and more as the axial contribution 
is enhanced. The S ww statistic measures the average of 
the modulus of the products of the sin of angles between 
quadrupole vectors, octopole vectors and the cos of the 
angle between the area vectors. Therefore on top of the 
information already contained in S nn the S ww statistic is 
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FIG. 6: Contour plot of the scatter of pairs (Sww , Sfu!i nc ) a fter 
a strong axial contingent of a to = 1000/iK is induced to the 
multipole vectors (see also Figures [1I2|I . The total number of 
Monte Carlo pairs is N = 10 J . The all-sky power statistic 
reacts heavily as the scale on the Sf^jJ nc -axis is shifted by four 
orders of magnitude with respect to the case of aeo = 70/xK 
(Figure[5|. The likelihood maximum departs very articulately 
from the WMAP(3yr) data. The contour lines are defined like 
in Figure [3] 



able to go to zero for highest zonal contamination as the 
closeness of the multipole vectors in that case dampens 
the product of sines and cosines quadratically to arbi- 
trary small values. Thus we find that S ww is the more 
convenient statistic for further analyses, as it does con- 
tain more information than the S nn statistic and addi- 
tionally shows a simple and clear asymptotic behavior. 
In the case of this statistic the anomaly is significant at 
99.5%C.L. with respect to S^ (3yr) = 0.7604. Similarly 
to before the maximal improvement is reached with an 
axial contribution of uio = 70^K, which degrades the 
anomaly in S ww to 99.2%C.L. 

Now we return to the correlation analysis of the align- 
ment with the pure multipole power Ci. When intro- 
ducing an axial effect, say ago — 100/xK, we improve the 
fit to the S ww statistic, but interestingly the multipole 
power anomaly becomes much more pronounced. This 
behaviour is expected [22, 23] for the C^-distribution (be- 
ing a modified ^-distribution) when the axial contribu- 
tion is enhanced, but it is unexpected that exactly the 
same happens for a multipole power distribution 'that 
knows of the intrinsic alignment of quadrupole and oc- 
topole'. This indicates that there is no correlation at all 
between multipole power and the phase alignment even 
when they are tuned to each other. 

Proceeding with the analysis of correlations between 
alignment and the full-sky power statistic, again we try 
to provoke correlation with the help of axial symmetry 
in the CMB. In fact we apply an axial effect of the ideal 
magnitude (a^o = 70//K) in order achieve larger values in 
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S W w The negative result is shown in Figure [5] As S'f™! 110 
is a linear combination of squared C# distributions it is a 
sharply peaked x 2 -hke distribution being very sensitive 
to axial contributions. Therefore the contour in Figure [5] 
is fairly shifted to the right (to higher values in Sf^Jf 10 ) 
and broadened with respect to the axially unmodified 
case, obviating any correlation with the intrinsic align- 
ment. The shape of the overall contour is roughly left 
invariant by the scale shift in S^ 110 . 

Figure [5] illustrates the pure zonal case. Here a whole 
ai — 1000/jK has been induced into the multipole vec- 
tors. Again, due to the sensitivity of S^jj 110 to axial con- 
tamination this pushes the allowed region in the scatter 
plot to very high values in full-sky power squared, de- 
generating the contour to a 'small' area far away from 
the measured three-year WMAP values. No change in 
correlation is observable. 

Obviously, no coupling of the multipole power statis- 
tic and the intrinsic alignment can be driven in favor of 
the anomalous experimental CMB data by an additional 
axisymmetric effect on top of the primordial fluctuations. 

V. CONCLUSIONS 

We have shown that a literal interpretation of the 'Axis 
of Evil' as an axisymmetric effect is highly incompati- 
ble with the observed microwave sky at the largest an- 
gular scales. The formalism of multipole vectors was 
used to separate directional information from the ab- 
solute power of multipoles on the CMB sky. Consid- 
ered were two choices of statistic, measuring the intrinsic 
cross- alignment between the quadrupole and octopole: 
the S nn and the S ww statistic. We confirm that the 
S ww statistic contains more information on the multi- 
poles and that it has more discriminative power as an 
axial effect is included. The presence of an axial symme- 
try in the CMB would excite zonal modes which are, in 
the frame of the axis, additional ago contributions in the 
language of the harmonic decomposition. Both statistics 
{S nn and S ww ) reach slightly better agreement with the 
measured values from the ILC(3yr) map at amplitudes of 
roughly ago = 100/iK. Further enhancement of the axial 
effect only reduces consistency with WMAP(3yr) data. 

Especially we have assayed in what way the alignment 
anomaly between quadrupole and octopole can affect the 
respective multipole power. We made several tests where 
we identified and selected the 'anomalous ag rn ' that are 
still consistent with data and checked whether the re- 
sulting distribution from these ag m for either power or 
alignment shows any change with respect to the unbiased 



case. For the all-sky multipole power we make use of the 
statistic S^™ - We demonstrated that the correlation 
between S^™ 13 and intrinsic alignment is only marginal 
(correlation coefficient of —0.0027). Thus a factorisation 
of the probability for the joint case into a product of the 
respective probabilities is allowed. 

We argued that the combined case of the measured 
all-sky power and the quadrupole-octopole alignment is 
anomalous at > 99.95% C.L. with respect to the WMAP 
three-year data. The correlation picture leaves no space 
for an axisymmetric effect in the large-angle CMB. 

These findings complement our previous studies [22j of 
the interplay of an axisymmetric effect and the extrinsic 
CMB anomalies (correlation with the motion and orien- 
tation of the Solar system Q). In that work it was shown 
that an axisymmetric effect might help to explain a So- 
lar system alignement. Finally, this study rules out that 
possibility. 

But there is a loophole. Here and in [22j we only con- 
sidered additive modifications of the ag m . Still, a pre- 
ferred axis could also induce multiplicative modifications 
in all ag m [l5| . This could avoid the problem of additional 
multipole power. However, multiplicative effects could 
only be achieved by non-linear physics, like systematics 
of the measurement or the map making process. 

A modelling that would be able to consistently remove 
both the power and the intrinsic alignment problem for 
low-£ must mobilize a more complex pattern of modifi- 
cations than the one induced by an axisymmetric effect. 
As already indicated by e.g. the odd extrinsic alignment 
with the ecliptic @, H, [|3, [ll, [26| the CMB anomalies do 
rather require a special plane than a preferred axis. The 
so called 'Axis of Evil' appears as just the normal vector 
of that plane, but no axial symmetry is present within 
that plane. 
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